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Abstract. Stellar oscillation under the combined influences of incident gravitational 
wave and radiation loss is studied in a simple toy model. The star is approximated 
as a uniform density ellipsoid in the Newtonian gravity including radiation damping 
through quadrupole formula. The time evolution of the oscillation is significantly 
controlled by the incident wave amplitude h, frequency v and damping time r. If a 
combination Hvt exceeds a threshold value, which depends on the resonance mode, 
the resonant growth is realized. 
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1. Introduction 

Stellar oscillation coupled to gravitational waves is one of important issues for the 
gravitational wave observatories such as LIGO, VERGO, GEO600 and TAMA. The 
current status of these detectors is given elsewhere(e.g, PQ). Detection of gravitational 
radiation from compact stars would provide valuable information about their interior. 
The frequency and damping time are uniquely affected by the stellar structure. 
Gravitational wave asteroseismology in the future is also discussed(e.g, 00] )• Various 
oscillation modes in neutron stars have been studied theoretically so far, but little is 
known of the excitation and damping processes. Most likely event that may excite the 
oscillations is the birth of the neutron stars after supernova. Some oscillation modes 
may also be induced by instabilities. An accurate description of these processes requires 
numerical solution of coupled system of hydrodynamics and the Einstein equations. 
Recently, the numerical simulations have been developing in the several approaches, 
e.g, 2D Newtonian simulation [4J, 3D smooth particle code simulation 0, Newtonian 
MHD simulation |6J, GR simulation |7] and the referees therein. 

Approximate approach is also helpful to get some insight into the physical processes. 
For example, excitation of neutron star oscillations in a close binary is studied by using 
following approximation [HI 13 EE IU|- One of the two stars is an extended body, and 
the other is approximated by a pointlike particle. The stellar perturbations in general 
relativity are solved for the extended body, while the pointlike mass moves on geodesic 
of the spacetime around the star. The oscillation mode can be excited by the particle, 
whose orbit is close enough to satisfy the resonance condition between normal mode 
frequency of the star and orbital Kepler frequency. Detailed discussion concerning 
excitation of neutron star oscillation in this approach may be found in Ref.[llJ and 
references therein. This is one of examples showing that oscillations are excited by 
the external almost periodic disturbance. The resonant growing mode is in general 
significantly affected by viscous damping, which may come from physical origins or 
numerical truncation errors. It is therefore better to examine resonant property in a 
simplified model. In this paper, we consider an idealized situation to study the resonant 
behavior. Incident gravitational waves are regarded as external periodic disturbance. 
The purpose of our study is the following. What is the resonance condition? How 
does the amplitude grow? The external disturbance is characterized by the frequency, 
duration and amplitude of the waves. These quantities are closely related with the 
resonance condition and the amplification of growing mode. The resonance condition is 
a relation between intrinsic and external frequencies. The amplification factor depends 
on growth rate of the mode and duration of the disturbances. 

This paper is organized as follows. In section 2, we briefly summarize the resonance 
appeared in a harmonics oscillator, which is helpful to examine the stellar model. In 
section 3, we consider the stellar oscillation with ellipsoidal model under Newtonian 
gravity. The gravitational radiation loss is included by the quadrupole formula. We 
show the numerical results in section 4. Section 5 is devoted to discussion. 
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2. Resonance in a harmonic oscillator 

We here briefly consider the resonance in a harmonic oscillator, which is a useful model 
to understand the more realistic case below. There are two kinds of resonance (e.g, 
[T2]). One type occurs when the frequency (= Q/2tt) of external perturbation matches 
with the intrinsic frequency (= u/2tt), i.e, Q = u. The relevant equation is given by 

X + u 2 X = /cos(fit). (1) 

In this case, the amplitude X increases linearly with the time, X oc (f/2Q)t sin(fit). 

Another type of resonance is parametric resonance, which occurs for time-dependent 
frequency. One example known as the Mathieu's equation is 

X + uo 2 (l + bcos(nt))X = 0. (2) 

The exponential growth for instability appears when uj/VL = n/2,(n = 1,2,3, •••). 
There is allowed range of frequency for each resonance condition, but the width of 
resonance range decreases rapidly with increasing n. The matching for large n is hardly 
satisfied. The strongest instability occurs for the fundamental one n — 1. Therefore the 
parametric resonance with n = 1 is very important. The growth rate s in X oc exp(st) 
for the fundamental mode can be calculated for small b as 



1 

S = -2 



(Q - 2uf + -b 2 uo 2 



1/2 1 

~ -boj (3) 



for the unstable region 



1 1 

boj < £1 - 2cu < -ho. (4) 

2 2 v ; 

In the presence of frictional damping, the amplitude decreases with time as exp(— jt). 
The damping counteracts the parametric resonance, so that the condition is more 
constrained. The growing solution of the resonance is possible only when b exceeds 
a threshold Aj/ui. 

3. Stellar pulsation 

Now we consider the pulsation of the star based on Newtonian gravity. The equation of 
motion for the pulsation driven by external acceleration g is 

dv 1 - - , 

— = —Vp-V& + g. 5 
at p 

As an external perturbation, we consider plane parallel wave with amplitude h and 
angular frequency Q, propagating along z-axis. The tidal acceleration by the incident 
gravitational wave is expressed as|13j 



9 



-hfl 2 x cos Q(t — z), -hQ 2 y cos Q(t — z), 



(6) 
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We here consider + mode only. Another polarization mode x is given simply by 7r/4 
rotation about z-axis, and may be ignored without loss of generality. Deceleration due 
to quadrupole radiation loss is given bv[T3] 

9 ^l:;\rV. (7) 



5c 5 " 



(5) 

where X\j is five time derivatives of reduced quadrupole moment. We incorporate both 
effects by eqs.flHJ and (JJJ). 

We are interested in global motion of a star, and therefore approximate it by the 
ellipsoidal model. In the approximation, the dynamical motion is limited to uniform 
expansion and compression along three axes. Thus, the fluid motion is described by 
three functions of time aj(i)(i = 1,2,3). The ellipsoidal approximation is often used 
as a simplified model to extract dynamical features in more realistic cases. (See e.g, 
|14 | IT5 | IT6 t lT7 t ll8j.) Assuming incompressible fluid with uniform density p, the dynamics 
can be determined by 

K 1 2G 

a 1 = - 2irGpa 1 A 1 H hQ 2 a 1 cos fit ^lf}ai, (8) 

a\ 2 5c & 



K 1 2G 
a 2 = - 2nGpa 2 A 2 + — + -hQ 2 a 2 cos fit - —2^a 2 , (9) 

a 2 2 bc° yy 

K 2G, 



03 = - 27rGpa 3 A 3 + - - -^Og, (10) 

a 3 5c° 



where the function Ai depends on the shape of the ellipsoid, and is described in pi 
The term K associated with pressure is algebraically determined from the condition 
d(\n(aia 2 a 3 )) / dt = 0. The system of eqs.(J5j)-(fTD|) depends on three dimensionless 
parameters, Q = ujt/2-k = ur, u/Q and h, where v[= u/2tt) is intrinsic frequency 
of the stellar oscillation, and r is damping time due to gravitational radiation. The 
amplitude becomes e -1 after Q times oscillations. These values are physically related 
with mass M and radius R of the star as 

/4GM\ 1/2 , 2GMR 2 uj 4 , , 

In order to clarify the dynamics of this system, we assume small deviation from 
spherically symmetric equilibrium (a, = 1), and linearize eqs.([H|)-([lL)p. We neglect the 
radiation loss for simplicity and limit the oscillation to the toroidal motion in x-y plane. 
The dynamics can be expressed by a single function X as 

h rfi\ 2 „ \ „ l 



X + u z U + 7r(jjJ cos fit J X = —-jhVt 2 cos fit, (12) 

where X represents the deviations 5ai{= a^ — 1) from the equilibrium state: 

8a\ = —5a 2 = X, 5a 3 = 0. (13) 

Equation (|T2*|) is the Mathieu's one with periodic source term. As considered in section 
2, it is clear that two kinds of resonance are possible in eq. ()12j) . and hence in the system 
of eqs. (jSjl-lfTUjl. In the astrophysical situation, periodic external perturbation may be 



Resonant growth of stellar oscillations 



5 



expected, e.g, in close binary system, where the standard resonance due to the periodic 
driving force is relevant. On the other hand, we do not definitely know the situation in 
which the parametric resonance occurs, but the relevant elements may be involved in 
the complicated system like the hydrodynamical core collapse. 

4. Numerical calculation 

4-1- Off-resonant case 

Before examining the resonant properties of the system (J%|-pO)l. we first of all show 
the behavior of non-resonant case. This is necessary for checking and comparing 
models. As for initial conditions of numerical calculations, we adopt static state, i.e, 
(ai = a2 = a>3 = 1, cl% = 0,2 = 63 = 0). The oscillation amplitudes Sa±, 5a% as a 
function of time are shown in Fig.l. The amplitude Sa2 is always well approximated as 
Sa2 ~ —Sa\, and is omitted in the figure. The parameters adopted for the calculation 
are Q = 90, uj/Vt = 0.6 and h = 0.1. The damping time r of the model corresponds to 
t = 2itQ/uj = 150 x (2ir/Q). In the early stage of the oscillation, the temporal profiles 
are a mixture of two modes with frequencies u(= 0.6Q) and Q. There is a modulation 
of two modes, which causes larger amplitudes of the oscillations due to interference 
of two waves. As the time goes on, the amplitudes of both modes are damped by 
the radiation loss. Since the external perturbation with Q is always supplied, the mode 
becomes dominant around r (= 150 x (2tt/Q)), and the oscillation is eventually enforced 
to match it at the late phase, typically at several times r, as shown in right panel of 
Fig.l. From the numerical results, we found that the motion along 03 is also induced 
by the incompressible condition as = 1/(0102). This differs from linear analysis, in 
which a 3 is always constant. The oscillation of 5a 3 is of order h 2 in the amplitude, and 
is induced by a mixture of stellar intrinsic oscillation, uo and external perturbation, Q. 
The frequencies of Sa 3 are therefore described by 2Q, 2u and Q ± uj. 




2 4 6 8 10 150 152 154 156 158 160 990 992 994 996 998 1000 

£2t/2n Qt/2n at/2lz 



Figure 1. Oscillation amplitudes Sa± and 60,3 as a function of time Ot/(27r). 
Parameters of the model are Q = 90, w/Cl = 0.6 and h = 0.1. From the left to 
right panels, the amplitudes of the oscillations are shown for the early, middle and late 
stages. 

We have also calculated the oscillations for different parameters, and found that 
the general behavior is almost the same unless the resonance condition is satisfied. 
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The temporal behavior may be easily understood by analyzing the numerical data with 
Fourier transformation. It is found that the oscillation is approximated by the following 
curve consisted of two frequencies uj and O: 

Vt 2 h 

Sai ~ — Sa 2 ~ —r^ TT (cos fit — exp(— t/r) cos cut) , 5a% oc h 2 . (14) 

2\IV — uj 2 \ 

This formula can be derived from the linearized system eq.(|12|) and the accuracy is 
checked by fitting the numerical results for a wide range of parameters. 

In the limit of t —>■ oo, the oscillation with uj in a\ or a-i is completely damped, and 
the mode matched with the incident wave survives with the amplitude Q 2 h/(2\Q 2 — uj 2 \). 
The amplitude of the steady state realized at t — > oo becomes larger when Q — > uj. The 
amplitude apparently diverges at the resonance, but does not actually, because even 
small dissipation is important and suppresses the divergence in that case. See section 
4.3. 

4-2. Parametric resonance of fundamental mode 

In this section, we consider the parametric resonance of the fundamental mode by setting 
uj/Q = 1/2. In Fig.2, we show time evolution of the model with h = 0.1, Q = 45. The 
damping time r corresponds to r = 90 x (2ir/Q). In the early phase of the time evolution, 
the amplitude of intrinsic stellar oscillation mode uj{= Q/2) exponentially grows up. 
After a hundreds of the oscillations, damping effect becomes important around t ~ r. 
Eventually the mode uj is damped, and oscillation is enforced to external mode with 
Q at the late phase t 3> r, in which the amplitude of 8a\ becomes 8a\ ~ 2h/3. This 
behavior in the final steady state is the same as that in off-resonant case. Compare 
the right panel in Fig.l (t w 1000 x (2vr/fi)) with that in Fig.2 (t w 700 x (2ic/ti)). 
Both figures show the behavior corresponding to several times r. Thus, the parametric 
excitation decays away in this case. 




2 4 6 8 10 90 92 94 96 98 100 690 692 694 696 698 700 
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Figure 2. Oscillation amplitudes Sai and 6123 as a function of time fit/(27r). 
Parameters of the model are Q = 45, u)/Cl — 0.5 and h = 0.1. 

We consider a case in which the radiation loss is less effective. The results for the 
parameters h = 0.1, Q = 75 are shown in Fig.3. The damping time r corresponds to 
t = 150 x (27r/0). In the early stage, a mixture of two frequencies uj(= Q/2) and Q can 
be seen. There are no drastic differences between the left panels in Fig.2 and Fig.3 as 
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for the early stages t < r. By comparing the right panels in Fig.2 and Fig.3, it is found 
that the parametric excitation occurs in Fig.3, since the oscillation amplitudes of both 
8ai and Sa 3 are indeed enhanced. The dominant oscillation frequency at the later stage 
is determined by not external one Q, but intrinsic one u. The frequency of £03 is also u, 
which comes from quadratic coupling between Q and u, as Q— u (= uj). The other higher 
frequencies such as 2u, VL + to (= 3u) and 2Q (= Au) are less excited. The amplitude 
at the later stage is saturated around the value 8a\ ~ 5a 3 « Ah, which is several times 
larger than that of off-resonant oscillation. This saturation property is quite different 
from the linear theory, in which the amplitude of the unstable mode exponentially 
increases so far as resonant condition is satisfied. In the nonlinear ellipsoidal model, 
the motion along the longitudinal direction 03 plays a key role on the saturation as 
explained below. The amplitudes of 5a\ and 5a2 grow in proportion to h, while that 
of Sa 3 , which is formally ~ h 2 , also grows due to nonlinear coupling. The amplitude 
of Sa 3 is given by 5a 3 ~ (5ai) 2 . When the oscillation amplitude becomes large enough, 
say 5a\ ~ 0.5, 5a 3 is also large enough. The motion of 8a 3 is no longer neglected and 
significantly affects the dynamics of 5a± and 5a 2 . In presence of the dynamics of £03 the 
growing amplitudes are saturated around the same order Sa± ~ Sa2 ~ Sa 3 ~ 0.5. 




2 4 6 8 10 150 152 154 156 158 160 1190 1192 1194 1196 1198 1200 
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Figure 3. Oscillation amplitudes 5a\ and 8a^ as a function of time fit/(27r). 
Parameters of the model are Q = 75, w/f2 = 0.5 and h = 0.1. 

So far we have shown two different results of the resonance in Fig.2 and Fig.3. The 
parametric oscillation is excited for weak damping case, i.e, large Q, while it decays 
out for strong damping case, i.e, small Q. The excitation condition should also depend 
on the amplitude h of the incident gravitational wave. In order to study the resonance 
condition for uj/Q = 1/2, we have performed the numerical calculations of eqs.(f%|)- ()10j) 
for various parameters (h,Q), typically in the range 20 < Q < 200. The results are 
summarized in Fig. 4. 

We found that a combination of the parameters, hQ is a very important factor to 
determine the evolution|. As shown in Fig. 4, there is a critical curve discriminating 
between growth and damped cases. The curve is empirically given by hQ ~ 7. Below 
the critical curve, the damping effect is so strong that the parametric excitation is 
suppressed like in Fig.2. On the other hand, above the critical curve, the damping 

f Note that hQ corresponds to frw/7 in the harmonics oscillator shown in section 2. The importance 
of this factor is reasonable. 
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effect is less effective. In this case, parametric excitation of the intrinsic oscillation with 
uj occurs and the amplitudes can grow up. Eventually, the nonlinear effect becomes 
important and the amplitudes of the unstable oscillations are saturated around finite 
values. 

0.2 
0.15 
h 0.1 
0.05 


20 40 60 80 100 120 140 160 
Q 

Figure 4. Excitation of the resonant oscillation significantly depends on the 
parameters, Q and h. Open circles and filled squares respectively denote growth and 
damped cases in the final steady states. The critical curve discriminating two states is 
approximately given by hQ « 7. 

4-3. Resonance at u = Q 

In this section, we consider the resonant behavior at uj = Q. Two typical cases are 
shown in Fig.5 and Fig.6. These figures show their long term evolutions up to hundreds 
of cycles of the oscillations, so that each sinusoidal curve is not clear, but the gradual 
change of the envelope is seen. The model parameters of Fig.5 are Q = 100, h = 5 x 10~ 3 , 
while those of Fig.6 are Q = 100, h = 5 x 10~ 2 . The behaviors are the same only for 
the initial stages, in which the amplitude of 5a 1 (or 5a 2 ) increases linearly with time: 
5ai pa hQt sin(Qt) / 4, so far as 8a\ <C 1. The later behaviors significantly depend on the 
adopted parameters, and are exhibited in quite different way. We have also calculated 
in a wide range of parameter space (Q, h) and found that the most important factor is 
hQ like in the previous section. The reason is explained below. 

The oscillation of 5a^ is inevitably induced by the nonlinear coupling: 5a 3 oc 
8a\5a2 ~ (hQt) 2 sin 2 (fit)/16. The amplitude of 5a% grows quadratically with time, 
and is therefore small only in the initial evolution. There is an epoch in which both 
amplitudes of 5a 1 and 5a 3 become the same order. This non-linear timescale t n is 
estimated as t n = 4/(M7) by setting 8a\ = 603, if the damping may be neglected. It is 
important to compare this timescale with that of the radiation loss r = Q x 2n/Q. If 
t < t n , which corresponds to hQ < 2/n pa 0.6, the radiation loss controls the evolution. 
As shown in Fig.5 (hQ = 0.5), the grown amplitude of 5a\ is saturated around the value 
5ai ~ hVtr ~ hQ. The amplitude of 5a^ grows up in much longer nonlinear timescale, 
which is not exactly but roughly given by t n . We have shown in Fig.5 the time evolution 
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for t < 250 x (27r/Q), but confirmed that the subsequent evolution for t > 250 x (2ir/Q) 
becomes almost steady state with some small modulation in the amplitudes. 

On the other hand, if r > t n , i.e, hQ > 2/ir, the nonlinear coupling is important. 
The behavior in this case is highly unstable as shown in Fig.6 (hQ = 5). The amplitudes 
5ai {i = 1,2,3) become 0(1). The nonlinear timescale in this model corresponds to 
t n ~ 13 x (27r/0), where the amplitude 5a 3 becomes large. The time evolution exhibits 
chaotic behavior^, and therefore the actuate numerical integration is rather difficult for 
much longer time, say, up to 10 3 x (2ir/Q) in Fig.6. 




0.8 - 



1.2 " 




0.8 " 

0.7 ' 1 1 1 1 

50 100 150 200 250 

Qt/2rc 

Figure 5. Oscillation amplitudes 5a\ and 8az as a function of time Vit/(2it). 
Parameters of the model are Q = 100, u = Cl and h = 5 x 10~ 3 

5. Summary and discussion 

In this paper, we have examined resonant oscillations in ellipsoidal stellar model with 
incoming gravitational wave. Resonance occurs when a certain matching condition holds 
between the stellar oscillation and the external periodic disturbance. In the presence of 
radiation loss, the evolution significantly depends on a combination hQ = hvT consisted 
of the wave amplitude h, frequency v and the damping time r. The parametric resonance 
of the fundamental mode can grow up only if the quantity hQ is larger than a critical 

§ In general, it may be necessary to measure e.g, Lyapunov exponent in order to judge whether the 
nonlinear behavior is chaotic or not. We did not calculate the exponent, but found that the numerical 
behavior is very sensitive like chaotic system. 




20 40 60 80 100 120 140 



at/2n 

Figure 6. Oscillation amplitudes Sa\ and Sa^ as a function of time Qt/(2ir). 
Parameters of the model are Q — 100, u = fi and h = 5 x 10~ 2 

value 7). The amplitude of stellar oscillation increases, but is saturated due to non- 
linear coupling, by which energy is transfered to the additional oscillatory motion. This 
point differs from linearized system, in which the amplitude exponentially grows up. If 
the quantity hQ is smaller than the critical value, the resonant oscillation is damped. 
The resonant behavior at uj = Q is also determined by the parameter hQ. If hQ is 
larger than a critical value 1), then the system exhibits chaotic behavior. We have 
also performed numerical calculations to explore the possibility of parametric resonance 
in higher modes uj/VL = n/2 (n = 3,4, • • •), and found that the higher modes are much 
more difficult to be excited. 

In our idealized model, the external perturbation is assumed to be constant. This 
corresponds to the situation that duration t c of the driving source is much larger than 
radiation damping time r, or non-linear time scale t n . If this condition does not hold, 
then the environmental effect would appear first in the evolution. For example, when 
the external perturbation stops at t c , the enhancement by the resonance also stops there, 
and the amplitude is subsequently damped around r(> t c ). 

We now consider the astrophysical relevance of the resonant oscillations. The f- 
mode oscillation of a neutron star is estimated as the frequency v = a few kHz, the 
damping time r = a few xlCT 1 s, and therefore Q 10 2 . In this case, the exotic 
case needs large amplitude of the gravitational wave h > Q^ 1 m 10~ 2 . The condition 
is realized only in hydrodynamical regime, e.g, in core collapse of supernova and in the 
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final phase of binary coalescence. The duration of periodic disturbance is also important 
factor. Typically, the timescale is of order 10~ 3 s, and may be too short for the mode 
to excited. Anyway, more realistic simulation is needed. 

Our result may be scaled to the excitation of the f-mode in white dwarfs. In this 
case, the frequency and damping time are given by v ~ 0.1Hz, r ~ 10 11 s, and therefore 
Q ~ 10 10 . In this case, the condition for the exotic case becomes less severe h > 1CT 10 . 
The region available for this large amplitude may be estimated by the total amount of 
the radiated energy AMc 2 . The distance from the periodic gravitational wave source is 
still limited to r ~ GAM/ (he 2 ) < lpc. 
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